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In our former study (J. Phys. A: Math. Theor. 43, (2010) 325210 or larXiv:1002. 0999H [math- 
ph]) we introduced a modified Fourier-Cattaneo law and derived a non-autonomous telegraph- type 
heat conduction equation which has desirable self-similar solution. Now we present a detailed in- 
depth analysis of this model and discuss additional analytic solutions for different parameters. The 
solutions have a very rich and interesting mathematical structure due to various special functions. 
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The heat equation propagates perturbation with infinite velocity, which is a well-known theoretical problem from a 
' long time. A historical review of the different way-outs can be found in [IHH, According to Gurtin and Pipkin 0-|3> 
the most general form of the flux in linear heat conduction and diffusion related to the flux q expressed in one space 
^f*) • dimension via an integral over the history of the temperature gradient 

where Q(t — t') is a positive, decreasing relaxation function that tends to zero as t — t' — > oo and T(x,t) is the 
: temperature distribution. 

There are two notable relaxation kernel functions: if Q\(s) — kS(s) where S(s) is a Dirac delta "function", then we 
will get back the original Fourier law. 

In out former study [|| we introduced the Q(t — t') — ( t _^J~ +LO y time dependent kernel which leads to a non- 
autonomous telegraph-equation 
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\^ . which has self-similar solution of the form of 
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T{x,t)=t-«j(^) :=t- a f(v). (3) 



The similarity exponents a and f3 are of primary physical importance since a represents the rate of decay of the 
magnitude T(x,t), while /3 is the rate of spread (or contraction if ft < ) of the space distribution as time goes on. 
k>( Substituting this into ©: 

f"(r 1 )t- a - 2 [ef3 2 r 1 2 ] + 
f{ri)rjt- a - 2 [eccp - e(3(-a - f3 - 1) - f3a] + 
f{r,)t- a - 2 [-ea(-a - 1) - aa] = /"(^)r ^ , (4) 

where prime denotes differentiation with respect to rj. 

One can see that this is an ordinary differential equation(ODE) if and only if a + 2 = a + 2/3 (the universality 
relation). So it has to be 

(3 = 1 (5) 

while a can be any number. The corresponding ODE is: 

f"(n)[e V 2 - 1] + f'(ri)r,(2ea + 2e-a) + f(n)a(ea + e-a) = Q. (6) 

In pure heat conduction-diffusion processes - no sources or sinks - the heat mass is conserved: the integral of T(x, t) 
with respect to x does not depend on time t. For T(x, t) this means 

J T(x, t)dx = t- a J f(j)dx = t- a+1 J f(r])dr) = const (7) 



-2- 



FIG. 1: Solutions for (Eq. 13) for zero and positive integer a/2e values. Thin solid line is for 0, thin dotted line is for 1, thick 
dashed line is for 2 and thick dash dotted line is for 3. 

if and only if a = 1. Eq. 6 can be written as 

(efv 2 /)" = afrf)' (8) 
which after integration and supposing f(rjo) = the first integration constant is zero for some 770 gives 

df arjdrj 



f erj 2 — 1 

The solution which is globally bounded and positive in the domain {(x,t) : 1 — erj 2 > 0} and has the form 



(9) 



f = (l-e V 2 )i' 1 (10) 

(With substitution and direct derivation it can be seen that / = (erj 2 ~ 1)^ _1 is a solution also.) The corresponding 
self-similar solution is 

T(x,t) = -(l-e^Y (11) 
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which was presented in our former study [3| with a detailed analysis. 

If the first integration constant is not zero (ci 7^ 0) than the following first ODE should be integrated 

e(r! 2 f( V ))' - ar,f(r,) = f( V ) + Cl (12) 
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FIG. 2: Solutions for (Eq. 13) for negative integer a/2e values. Thin solid line is for -1, thin dotted line is for -2, thick dashed 
line is for -3 and thick dash dotted line is for -4 



after some algebra the solution can be written in the following form 
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Note, that the former solution is still there independently and as a kind of form factor. Where 2-F1 is the hypergeometric 
function [5 8]. It is well known that many elementary functions can be expressed via the hypergeometric function. 
Fortunately, in our case if ~ is an integer or a half-integer the hypergoemtric series breaks down to a finite sum of 
terms. There are four basic cases which has crucial importance: 
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for half-integer values 



a _ 1 
Ye ~ Y 
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With the following recursion relation all the other cases can be evaluated 

(c — a)2-F\(a — 1, 6; c; z) + (2a — c — az + 2)2^1(0, b; c; z) + a(z — 1)2^1(0 + 1, b; c; z) = 0. 
Even for negative parameters there are closed relations 
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FIG. 3: Solutions for (Eq. 13) positive half-integer a/2e values. Thin solid line is for 1/2, thin dotted line is for 3/2, thick 
dashed line is for 5/2 and thick dash dotted line is for 7/2. 
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Let's consider ci = 1 and c-i = solutions and analyze it in details. We will examine 4 different cases, (from 
physical reasons the propagation velocity e always have to be positive, we take the plus unity value, and changing the 
'a' parameter). For the first case let's take a/2e > with integer values. The resulting curves are presented on Figure 
1. We can see, that for positive integer a/2e values the solution is compact but have a finite jump at ±1. Figure 2 
shows the same results but when a/2e is equal to -1,-2,-3 and -4. All figures are continuous at the [-1:1] and have 
singular asymptotic values at the boundaries. Figure 3 presents the results for positive half-integer a/2e values. All 
the curves are odd functions. Figure 4 presents curves for negative half integer a/2e values. All the results are odd 
functions and have asymptotes at ±1. 

Figure 5 presents the self-similar solution 



T(x,t) 



1 



{signum(e(x/t) 2 - 1) & ){-signum{e(x/tf - 1))£) • (|) - 2 Fi Q, ^; ^;e{x/t) 2 
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for a/2e = 2. We can see the discontinuous jump at the x=t line. For large x and t the solutions goes to zero. 



(21) 



In the following we present solutions where alpha is not unity, we consider the a = 



—2 and /? = +1 values. (For 
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FIG. 5: The self-similar T(x,t) solution (Eq. [2T) for a/2e = 2. 



6 



arbitrary alphas the solutions are similar to this case but the details are much more complicated to understand.) Now 
(Eq. 4) becomes 

f'mer, 2 - 1] - f'( V ) V [2e + a] + 2f(rj)[e + a] = 0. (22) 
For general values of e and a the solutions are the following: 

f(r)) = c 1 P|+ 2 (^)(^ 2 - +c 2 Q j H(V~eri)(eV 2 ~ (23) 

where ^f^-ilv^ 7 ?) an d Qa/^e-iXv^ 7 ?) are the associated Legendre functions of the first and second kind [5l4l0|. 
This result can be obtained by tedious calculations transforming to the standard form of the Legendre differential 
equation Kamke [ll[ or by using Maple 9.1 program package. The first and second arguments of P and Q are the 
order v — a/2e — 1 and the degree /i = a/2e + 2 of the Legendre function, which may take unrestricted real values. 
The variable y = ^/erj may have complex values es well. 

It is well known that if the order and the degree are integer numbers (y = leN, [i = meN) and the order is 
larger that the degree (v > /i) the associated Legendre functions of the first kind became the associated Legendre 
polynomials. These polynomials span a Hilbert space with the following orthonormalization condition: 

- 2 h- W (24) 



2h + l(h- m)\ 



Similar relations are available for the second order associated Legendre polynomials Q™(x) as well. 

If the order and the degree are non-integer values than the Legendre functions can be evaluated with the hyperbolic 
functions in our case the formulas are the following: 

P ^ {V * v) (^-i)- +1 r(-f -i) 2 1 W ~^ ' " ' 2 " ~T 

and 



1 i(4e+ a )n / a \ / a 3a a 11 y/erj 



X 



2e 2' 2 2 

2(*+ 1 )(^ry-l)^- 1 (^+l)^V[e(^) f+2 ] ■ (25) 

The function 2-Fi(--) i s the hypergeometric function again, and T(...) is the gamma function [5, 6]. In our case, both 

the hypergeometric functions has a large degree of symmetry hence, the series expansion breaks down after some 
terms 

p (a a a 1 y/eiA . (1 + ^ + g^fKj + ffi ^ 

Ur> + 1, -tt: - i, o ~ - 4 /i , o /r„ i , Ia (26) 



2e' 2e ' 2e '2 2 j (1 + 2^lr] + erf){2e + a) 



and 



/ a 3 a a 1 1 ^77 \ e(l + e?? 2 + a?? 2 ) 

V2e 2 2e 2e 2 2 2 J (e + a)(e ^ 2 _ 1} - +1 

Now, after some algebra the total regular solution can be obtained as follows (ci = 1, c% = 0) 

M = Pt + t{yT*,)tf - 1)- +1 = ^+^ + W(^ (28) 
JK " -,- iyv " y 11 r (=f - 1) (2e + a) V ; 



Analyze now the results, which means four different cases again. (a/2e is positive/negative integer, or a/2e posi- 
tive/negative half-integer). 
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If 'a' is positive or negative odd number, then a/2e is a half integer number and the domain of (y/erj + i) a / 2t + 2 is 
for 77 > —1 which means that the (1 + r/ 2 (e + a)) parable is only considered for 77 > —1. Such solutions can be seen 
on Fig 6. 

If 'a' is an even number then the cases are a bit more difficult. From the properties of the Gamma function is clear 
that if a/2e is a negative integer, (a is a positive even number) than f(rj) is zero. So 'a' have to be negative and even. 
From the weight factor of (Eq. 23) (erj 2 — we can see that if a is divisible with 4 than the power of (e?/ 2 — 1) is 

an integer, and the range is the whole real axis, however if a is negative even but not divisible with 4 than the power 
of (erj 2 — 1) is a half-integer, which means the range is not continuous on the whole real axis. This property together 
with the Legendre function gives a cut at 77 = — 1. Such solutions are presented on Fig 7. 

The total self-similar solution (Eq. 3) is presented for a = —2, e = 1 




T(x,t) = ?f(x/t)=fPll-) K - -1 (29) 



of Figure 8. The solution is a nice continuous function. 

For the irregular solution (ci = 0, c-2 — 1) the solution formula is however a bit more complex, however still 
four different cases have to be analyzed [(a/2e) is positive/negative integer or positive/negative half-integer value]. 
The irregular Q Legendre function is not defined for any negative (integer and half-integer values also) order values 

(a/2e — 1). The complex exponential factor in (Eq. 25) e ( at can induce pure real, pure imaginary or mixed 
solutions. E.g. for ^ + 1 = 3/2 or 7/2 the next formula is valid 

Z(£77 2 - 1)£ + 1 = (l-£7 ? 2 )^+ 1 (30) 

which interchanges upper or lower boarders. 

Figure 9 presents the results for e = 1, a = 1, 2,3, 4. The domain of the (rj 2 — 1) 5 ^ 4 Q 2 _ 1 (77) is for 77 < —1 the range is 
also bounded from below. The second solution (x 2 — 1) 3 / 2 Qq(ie) has a continuous range and domain with a jump at 

77 = — 1. Solution (?7 2 — l) 7 / 4 <2?(77) is very similar to the first solution. The domain of the last curve (?/ 2 — 1) 2 Q\(ji) 

2 

is the — 1 < 77 < 1 opened interval and the function has singular values at the boarders. 




FIG. 8: The self similar-solution for a = —2, e = 1 
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FIG. 9: The results for the irregular Legendre function the solid line if for a = 1 and the dotted line is for a — 2 the dashed 
line is for a — 3 and the dot-dashed is for a — 4. 



At this point, we mention than all kind of solutions can be written in the form of the product of two travelling waves 
propagating in opposite directions. If we insert c 2 = 1/e (the wave-propagation speed) into an irregular solution e.g. 



T(x,t) = t 2 Q\ 



'ex 



(!)' 



1 



after some algebraic manipulation we get 

T(x,t) 



i 2 Qt(— ) U(x-ct)U(x + ct) 



(31) 



(32) 



where U(x ± ct) = (x ± ct) 2 . Which means a distorted wave solution with a non-trivial weight function. Figure 10 
presents a self-similar irregular solution (Eq. l3"Tj) . Note, the cut at the x = t line. We think that after this kind of 
analysis all the properties of the solutions are examined and understood. 



Let's consider now the 2 dimensional telegraph- type equation of the following form: 

c d 2 S(x,y,t) + a dS(x,y 7 t) = d 2 S(x,y,t) + d 2 S(x,y,t) 
dt 2 t dt dx 2 dy 2 

To avoid further confusion we are looking for the solution in the form of 



(33) 



*(■'■■ //./) = /■■'../ (^) :=/-"//('/')■ (■■'>4! 



We may consider this as a the L 1 vector norm as well. 

The idea behind this Ansatz is the exchange symmetry over the spatial coordinates: 



^feM = d'Sjxft) = t ^ w -v (35) 



dx 2 dy 2 
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where prime denotes differentiation with respect to r\. 

If we use the condition that (3 = +1 we get the most general ordinary differential equation: 

S"(?7)[«7 2 - 2] - g'(r))r]{a - 3e) + g(v)a(ae + e-a) = 0. 

The solution is the following: 

g(ri) = Cl (eri 2 -2)£-ip*~' 

c 2 {eif -2)£-tQZ 
where P and Q are the associated Legendre functions, and 



(36) 



2 V 



2" 



(37) 



Ji-ka. + 4 - 4a 2 )e 2 + 4a(a - l)e + a 2 - e , . 

w = (38) 

If a = +1 than |a| > |2e| is a constrain for a real u> number which is the order of the Legendre function. If we consider 
the —4a + 4 — 4a 2 = condition than it comes out a.\ — —y/E/2 - 1/2 or «2 = v5/2 — 1/2. This means for a>i 
the a < and for «2 the a > conditions have to be fulfilled to have a real order for the Legendre function. It is 
interesting to mention that the number ~ 1.61803. is the golden ratio. 

If j3 = +2 the ODE and the solutions are the followings: 

g"(v)[ev 2 - 2] - g'(vHa - 3e) + 2g(r,)(e + a) = Q. (39) 
g( V ) = Cl (er, 2 - 2)^ 1 /4pa/2A-i/2 ( ^/2 ?7) + Cl(e?7 2 _ 2)*- 1 /4Q»/a/e-i/a (v ^ ?) (40 ) 



where o has a bit simpler form: 
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to have real number for the order of the associated Legendre function the " a" and e values are not independent from 
each other, and the a/2e relation became a bit more complicated. However, the structure of the solutions remains 
the same, like for the 1 dimensional case. Some solutions are continuous on the whole real axis and some of them 
have a bonded range from above or from below. 

As for the more usual Euclidean norm we may consider the following Ansatz as well 



S(x, y, t) = t- a g y X2 t e V j ■= t~ a g(v) (42) 

where \/ ' x 2 + y 2 is the usual distance between the two points, which is never negative. After having done the usual 
derivation we got (3 = 1 and the following ODE 

g"(n)(en 2 - 1) + g'{rj){[2ea + 2e - a]n -l/rj) + g{n){ea 2 + ea - aa) = (43) 

the most general solution is the following: 

fa (l + q)e-a 2ae + e - a 2 
9(V) = ci -2 F\ ; ; 1 - er, 

c 2 - 2 F l \l--, ;i ^ ;1-^ 2 J (ery 2 - 1) * . (44) 

Which means that even this Ansatz is not contradictory, and gives a reasonable ODE. Note, that the only basic 
difference between this ODE and the original one (Eq. 4) is the extra 1/rj term after the first derivative of girf). From 
the (en 2 — a) form factor of the second derivative we can see, [3] that this ODE also has singular points at n = ±1 for 
e = 1 which is not special for wave-equations. If these singularities are relevant or can be eliminated depends on the 
relations of the parameters, (see Fig. 9. and explanations there). Here, we can find large number of analytic solutions 
again (for integer/half integer a, e and a), which are similar to (Eq. 14 - 20). 

We hope that with the help of these two mentioned models we can analyze 2 dimensional dissipative flows int the 
future where the the original telegraph equation plays an important role [12j . 



At last we investigate the one-dimensional telegraph-type equation with a general source term 

d 2 U{x,t) adU(x,t) d 2 U(x,t) 



e- 



dt 2 t dt dx 2 
We are looking for solution of (|4"5)l of the form 



+ s(x,t)U(x,t) =0 (45) 



U(x,t)=t~ a h[^) := t- a h(r,). (46) 

(To avoid confusion we use different letters for the basic variable, and for the Ansatz as well.) After the standard 
derivation and parameter investigation it comes out that j3 = 1 is a necessary condition. The most general ODE 
without a source term is still (Eq. 4.). Note, that any kind of source term (without derivation) gives additional terms 
for the last parenthesis which should be a function of r\. There are two important source terms in physics, the first is 
the harmonic oscillator. So s(x,t) — Dx 2 where D is the stiffness of the oscillator. This automatically fix a to +2. 
Our ODE is now: 

ti'(n)[er) 2 - 1] + h'(rf)ri[6e -a]- h(r})[-6e + 2a + Dr/ 2 } = 0. (47) 

the solution is the 

1 -4e + a D 10c - 3a 2 \ r 2 



%) = d ■ HeunC I 0, --, — ^— , - — , ^ , of ) ■ [erf - lj 2 « + 

„ ( n 1 — 4e + a D lOe - 3a , > , , . _ 
c 2 • HeunC ( 0, -, - ,- — , , tif ) V [en 2 - l] 2 < (48) 
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HeunC function which has three singularities. The regular ones are in T^tj and the irregular is in infinity. There 
are different kind of Heun functions, and all kind can be evaluated with a generalized power series, or for special 
arguments represents different kind of special functions like the hypergeometric functions, Legendre functions, Bessel 
functions, Gegenbauer polinoms, Kummer functions and so on. Unfortunately, for the above given parameters, the 
Heun function does not represent any kind of special functions, but the series expansion is still valid. 

_V 1 -4e + a D lOe - 3a 2 \ 1 / a\ 2 

ae a 2 c\ A ( a 3 (4e 2 + 2c)a ec\ fi ^. s . , . 



12 24 12/ ' V 240 240 20, 

The other important solutions is the s(x,t) — q/x case, which is the Coulomb potential where q is the electrical 
charge . This automatically hx a — — 1. Now the ODE is the following 

h"(n)[erf - 1] - ti(n) V a - h(n)[-a + q/rj] = 0. (50) 

The solution is 



„ ( ae + a\/e a a _ 

fcfa) = Cl • HeunG [2, ^ , -1, --, --,0, ^fen + 1 

c 2 h 2 -!]* [-Viry-1]^ [-v^ + l]"*- 
HeunG I 2, 2 A_ _ i _ 2 + -, 0, V^V + 1 J (51) 

This solution has four singular points at at oo and at ^fl/^/e. Without any physical restriction between the three 
parameters, further pure mathematical investigation has not much sense. 

We think that this model can help us to investigate further the features of the quantum telegraph equation fl3l ], 
there the first time derivative has an additional complex unit parameter. 

In summary: 

In the presented study we gave an in-depth analysis of a time-dependent telegraph-type equation for heat propagation. 
All the cases for different parameter ranges are carefully examined and analyzed. Large number of special functions 
came into play, which may serve as a strong hint for further relevance of this equation. We hope that our equation 
can help us to investigate two dimensional turbulent flows, or even some quantum mechanical problems. 
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